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The detailed spectral and spatial characteristics of the signature imprinted on the 
cosmic microwave background (CMB) radiation by Compton scattering of the ra- 
diation by electrons in the hot gas in clusters of galaxies - the Sunyaev-Zeldovich 
(S-Z) effect - are of great astrophysical and cosmological significance. In recent 
years observations of the effect have improved tremendously; high signal-to- noise 
images of the effect (at low microwave frequencies) can now be obtained by interfer- 
ometric arrays. In the near future, high frequency measurements of the effect will 
be made with ground based and balloon-borne telescopes equipped with bolome- 
teric arrays. Towards the end of the decade the PLANCK satellite will carry out 
an extensive S-Z survey over a wide frequency range. Along with the improved 
observational capabilities, the theoretical description of the effect, and its use as 
a precise cosmological probe, have been considerably advanced. In this review, I 
briefly discuss the nature and significance of the effect, its exact theoretical de- 
scription, the current observational status, and prospects for the near future. 

1 Introduction 

The S-Z effect was first described by Zeldovich & Sunyaev (1969) and Sunyaev 
& Zeldovich (1972), and after many attempts over a decade, it was convincingly 
detected in several clusters by single-dish radio telescopes (for general reviews, see 
Rephaeli 1995a, Birkinshaw 1999). Increased realization of the cosmological signif- 
icance of the effect has led to major improvements in observational techniques and 
to extensive theoretical investigations of many of its facets. The most important 
recent development is the ability to carry out sensitive radio interferometric mea- 
surements of the effect (Jones 1993, Carlstrom et al. 1996). Impressive images of the 
effect in more than thirty clusters have already been obtained using the OVRO and 
BIMA arrays (Carlstrom et al. 1999, Carlstrom et aL2000). Theoretical treatment 
of the S-Z effect has also improved, starting with the work of Rephaeli (1995b), who 
performed an exact relativistic calculation and demonstrated the need for such a 
more accurate description. 

The effect is a direct probe of clusters; high resolution measurements yield 
the spatial distributions of the hot intracluster (IC) gas and the total mass, as 
well as information on the evolution of clusters. Of more basic importance is the 
ability to determine the Hubble constant, Ho, and the density parameter, f2, from 
S-Z and X-ray measurements. This method to determine Ho, which has clear 
advantages over the traditional galactic distance ladder method, has already been 
employed but its full potential has not yet been realized because of substantial 
observational uncertainties. Sensitive mapping of the effect, and the understanding 
and minimization of systematic uncertainties in the S-Z and X-ray measurements, 
are essential steps towards this goal and constitute the main challenge of current 
and near future work on the S-Z effect. 
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The recent progress in S-Z work, and especially the improved spectral and spatial 
capabilities that will be attained in the near future by ground based and strato- 
spheric projects, will advance this field to the forefront of cosmological research. 
Here I review recent progress in this field and briefly discuss the prospects for the 
near future. 

2 Exact Description of the Effect 

Compton scattering of the CMB by hot gas heats the radiation, resulting in a 
systematic transfer of photons from the Rayleigh- Jeans (R-J) to the Wien side of 
the (Planckian) spectrum. An accurate description of the interaction of the ra- 
diation with a hot electron gas necessitates the calculation of the exact frequency 
re-distribution function in the context of a relativistic formulation. The calculations 
of Sunyaev & Zeldovich (1972) are based on a solution to the nonrelativistic Kom- 
paneets (1957) equation. The result of their treatment is a simple expression for the 
intensity change resulting from the scattering off electrons with thermal velocity dis- 
tribution, Alt, in terms of the Comptonization parameter, y = J (kT e /mc 2 )naTdl. 
This is a line of sight integral (through the cluster) over the electron density (n) and 
temperature (T e ); er<r is the Thomson cross section. AI t is negative in the R-J region 
and positive at frequencies above a critical value ~ 217 GHz (x = hv/kT = 3.83). 
Typically, in a rich cluster, y ~ 1CT 4 , along a line of sight through the center, 
and the magnitude of the relative temperature change due to the thermal effect is 
AT t /T = -2y 2 • 1CT 4 in the R-J region. 

The effect has a second component when the cluster has a finite (peculiar) 
velocity in the CMB frame. This kinematic Dopplcr component, which is obviously 
proportional to the component of the cluster peculiar velocity along the line of 
sight, v r , and to the Thomson optical depth, is Alk = x 4 e x (v r / c)t / (e x — l) 2 . The 
related temperature change is AT^/To = — (iv/c)t (Sunyaev & Zeldovich 1980). 

The nonrelativistic treatment of Sunyaev & Zeldovich (1972) is generally valid at 
low gas temperatures and at low frequencies. Rephaeli (1995b) has shown that this 
approximation is insufficiently accurate for use of the effect as a precise cosmologi- 
cal probe: Electron velocities in the IC gas are high, and the relative photon energy 
change in the scattering is appreciable enough to require a relativistic calculation. 
Using the exact probability distribution in Compton scattering, and the relativisti- 
cally correct form of the electron Maxwellian velocity distribution, Rephaeli (1995b) 
calculated AI t in the limit of small r, keeping terms linear in r (single scattering) . 
Results of this semi-analytic calculation, shown in Figure 1, demonstrate that the 
relativistic spectral distribution of the intensity change is quite different from that 
derived by Sunyaev & Zeldovich (1972). Deviations from their expression increase 
with T e and can be quite substantial. These are especially large near the crossover 
frequency (where AT t changes sign) which shifts to higher values with increasing 
gas temperature. 

The work of Rephaeli (1995b) sparked considerable interest which led to var- 
ious generalizations and extensions of the relativistic treatment of the S-Z effect. 
Challinor & Lasenby (1998) generalized the Kompaneets equation and obtained 
analytic approximations to its solution for the change of the photon occupation 
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Figure 1. The spectral distribution of AI t /r [in units of (he) 2 /(2(kT ) 3 )]. The pairs of thick (red) 
and thin (blue) lines show the relativistic and nonrclativistic distributions, respectively. Three 
pairs of lines are shown corresponding to kT e = 5, 10, and 15 keV. 



number by means of a power series in 9 e = kT e /mc 2 . Itoh et al. (1998) adopted 
this approach and improved the accuracy of the analytic approximation by expand- 
ing to fifth order in 9 e . Sazonov & Sunyaev (1998) and Nozawa et al. (1998) have 
extended the relativistic treatment also to the kinematic component obtaining - 
for the first time - the leading cross terms in the expression for the total intensity 
change (Al = AI t + Alk) which depend on both T e and v r . An improved analytic 
fit to the results of numerical integration (of the collision term in the Boltzmann 
equation), valid for 0.02 < 9 e < 0.05, and x < 20 (y < 1130 GHz), has recently 
been given by Nozawa et al. (2000). In view of the possibility that in some rich 
clusters r ~ 0.02 — 0.03, the approximate analytic expansion to fifth order in 9 e 
necessitates self-consistent inclusion also of multiple scatterings, of order r 2 . This 
has been accomplished by Itoh et al. (2000), and Shimon & Rephaeli (2001). 

The relativistically correct calculation has to be used in all high frequency S-Z 
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work, especially when measurements of the effect arc used to determine precise 
values of the cosmological parameters. Also, since the ability to determine pecu- 
liar velocities of clusters depends very much on measurements at or very near the 
crossover frequency, its exact value has to be known. This necessitates knowledge 
of the gas temperature since in the exact relativistic treatment the crossover fre- 
quency is no longer independent of the gas temperature, and is approximately given 
by ~ 217[1 + 1.167fcT e /mc 2 - 0.853(fcT e /mc 2 ) 2 ] GHz (Nozawa et al. 1998a). 

Use of the S-Z effect as a cosmological probe necessitates also X-ray mea- 
surements to determine (at least) the gas temperature. Therefore, a relativisti- 
cally correct expression for the (spectral) bremsstrahlung emissivity must be used 
(Rephaeli & Yankovitch 1997). In the latter paper first order relativistic correc- 
tions to the velocity distribution, and electron-electron bremsstrahlung, were taken 
into account in correcting values of Hq that were previously derived using the 
nonrelativistic expression for the emissivity (see also Hughes & Birkinshaw 1998). 
Nozawa et al. (1998b) have performed a more exact calculation of the relativis- 
tic bremsstrahlung Gaunt factor. 

Compton scattering of the CMB in clusters also affects its polarization towards 
the cluster. Net polarization is induced due to the quadrupolc component in the 
spatial distribution of the radiation, and when the cluster peculiar velocity has a 
component transverse to the line of sight, v± (Sunyaev & Zeldovich 1980, Sazonov 
& Sunyaev 1999). The leading contributions to the latter, kinematically- induced 
polarization, are proportional to (v±/c)t 2 and (v±/c) 2 t. Itoh et al. (2000) have 
included relativistic corrections in the expression for the kinematic polarization. 

3 Recent Measurements 

The dramatic improvements in the quality of S-Z measurements over the last seven 
years is largely due to the use of interferometric arrays. Dish arrays have several 
major advantages over a single dish, including insensitivity to changes in the at- 
mospheric emission, sensitivity to specific angular scales and to signals which are 
correlated between the array elements, and the high angular resolution that enables 
nearly optimal subtraction of signals from point sources. With the improved sen- 
sitivity of radio receivers it became advantageous to use interferometric arrays for 
S-Z imaging measurements (starting with the use of the Ryle telescope by Jones et 
al. 1993). Current state-of-the-art work is done with the BIMA and OVRO arrays; 
images of some 35 moderately distant clusters (in the redshift range 0.17 — 0.89) 
have already been obtained at ~ 30 GHz (Carlstrom et al. 1999, 2000). A beautiful 
example is the 28 GHz image of the cluster CL0016+16 obtained with the BIMA ar- 
ray (Carlstrom et al. 1999) which is shown in Figure 2. The ROSAT X-ray image is 
superposed on the contour lines of the S-Z profile in the lower frame; note the good 
agreement between the orientation of the X-ray and S-Z brightness distributions. 

Work has begun recently with the CBI, a new radio (26-36 GHz) interferometric 
array (in the Chilean Andes) of small (0.9m) dishes on a platform with baselines in 
the lm to 6m range. The spatial resolution of the CBI is in the 3' — 10' range, so 
(unlike the BIMA and OVRO arrays) it is suitable for S-Z measurements of nearby 
clusters. Some 9 clusters have already been observed with the CBI (Udomprasert 
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Figure 2. S-Z and X-ray images of the cluster CL0016+16. The S-Z image (false color) 
of the cluster, obtained with the BIMA array, is shown in the upper frame. In the lower 
frame, contours of the S-Z effect in the cluster are superposed on the ROSAT X-ray (false 
color) image (from Carlstrom et al. 1999). 
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et al. 2000). 

Recent work at higher frequencies includes measurements with the SuZIE ar- 
ray, and the PRONAOS and MITO telescopes. SuZIE (which was previously used 
to measure the effect towards the moderately distant clusters A1689 & A2163 
Holzapfel et al. 1997a, 1997b) observed the effect towards A1835 at three spectral 
bands centered on 145, 221, 279 GHz (Mauskopf et al. 2000). Observations at four 
wide spectral bands (in the overall range of 285-1765 GHz) were made of A2163 
with the PRONAOS atmospheric 2m telescope (Lamarre et al. 1998), leading to 
what seems to be the first detection of the effect by a balloon-borne experiment. 
The new 2.6m MITO telescope - which currently operates at four spectral bands 
and has a large ~ 17' beam - was used to observe the effect in the Coma cluster 
(Lamagna et al. , these proceedings). The strongest S-Z effect, with y = 1.2 x 10~ 3 , 
seems to have been detected towards the distant (z=0.45) cluster RXJ 1347 (Pointe- 
couteau et al. 1999) with the Diabolo bolometer operating at the IRAM 30m radio 
telescope. The Diabolo has a 0.5' beam, and a dual channel bolometer (centered 
on 2.1 and 1.2 mm). Four other clusters were also observed with diabolo (Desert 
et al. 1998). 

4 The Effect as a Cosmological Probe 

The S-Z effect is a unique probe whose indispensable qualities are being increasingly 
exploited for the determination of cluster and global cosmological parameters. For 
an extensive discussion of use of the effect as a probe, see the reviews by Rephaeli 
(1995a), and Birkinshaw (1999). 

In principle, high spatial resolution measurements of the S-Z effect yield the gas 
temperature and density profiles across the cluster. The deduction of an approx- 
imate 3D density and temperature distributions from their sky-projected profiles 
requires use of a deprojection algorithm. One such was developed by Zaroubi 
et al. (1998; see also the contribution by Zaroubi to these proceedings). Cluster 
gas density and temperature profiles have been mostly deduced from X-ray mea- 
surements. S-Z measurements can more fully determine these profiles due to the 
linear dependence of AI t on n, as compared to the n 2 dependence of the (thermal 
bremsstrahlung) X-ray brightness profile. 

The cluster full mass profile, M(r), can be derived directly from the gas density 
and temperature distributions by solving the equation of hydrostatic equilibrium 
(assuming, of course, the gas has reached such a state in the underlying gravitational 
potential). This method has already been employed in many analyses using X- 
ray deduced gas parameters (e.g., Fabricant et al. 1980). Grego et al. (2001) have 
recently used this method to determine total masses and gas mass fractions of 
18 clusters based largely on the results of their interferometric S-Z measurements. 
Isothermal gas with the familiar density profile, (l+r 2 /r 2 )~ 3/3 / 2 , was assumed. The 
core radius, r c , and (3 were determined from analysis of the S-Z data, whereas the 
X-ray value of the temperature was adopted. Mean values of the gas mass fraction 
were found to be in the range (0.06 — 0.09)ft. -1 (where h is the value of H in units 
of 100 km s _1 Mpc -1 ) for the currently popular open and flat, A-dominated, CDM 
models. 
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Measurement of the kinematic S-Z effect yields the line of sight component of 
the cluster peculiar velocity (v r ). This is observationally feasible only in a narrow 
spectral band near the critical frequency, where the thermal effect vanishes while the 
kinematic effect - which is usually swamped by the much larger thermal component 
- is maximal (Rephaeli & Lahav 1991). SuZIE is the first experiment with a 
spectral band centered on the crossover frequency. Measurements of the clusters 
A1689 and A2163 (Holzapfel et al. 1997b) and A1835 (Mauskopf et al. 2000) yielded 
substantially uncertain results for v r (170^, 490+^°, and 500 ± 1000 km s"\ 
respectively). 

Perhaps the most important use of the S-Z effect so far has been the mea- 
surement of the Hubble constant, Ho, and the cosmological density parameter, Q. 
Briefly, the method is based on determining the angular diameter distance, cIa, 
from a combination of Alt, the X-ray surface brightness, and their spatial profiles. 
Averaging over the first eight determinations of H (from S-Z and X-ray measure- 
ments of seven clusters) yielded Ho ~ 58 ± 6 km s _1 Mpc -1 (Rephaeli 1995a), but 
the database was very non-uniform and the errors did not include systematic un- 
certainties. Repeating this with a somewhat updated data set, Birkinshaw (1999) 
deduced essentially a similar mean value. A much larger S-Z data set is now avail- 
able from the interferometric BIMA and OVRO observations, and since the redshift 
range of the clusters in the sample is substantial, the dependence on f2 is appre- 
ciable. A fit to 33 cluster distances gives Ho — 60 km s _1 Mpc -1 for il = 0.3, 
and Ho = 58 km s _1 Mpc -1 for Q = 1, with direct observational errors of ±5% 
(Carlstrom et al. 2000). The main known sources of systematic uncertainties (see 
discussions in Rephaeli 1995a, and Birkinshaw 1999) are presumed to introduce an 
additional error of ~ 30% (Carlstrom et al. 2000). The current number of clusters 
with S-Z determined distances is sufficiently large that a plot of oIa vs. redshift (a 
Hubble diagram) is now quite of interest, but with the present level of uncertainties 
the limits on the value of Q are not very meaningful (as can be seen from figure 11 
of Carlstrom et al. 2000). 

CMB anisotropy induced by the S-Z effect (Sunyaev 1977, Rephaeli 1981) is 
the main source of secondary anisotropy on angular scales of few arcminutes. Be- 
cause of this, and the great interest in this range of angular scales - multipolcs 
(in the representation of the CMB temperature in terms of spherical harmonics) 
i > 1000 - the S-Z anisotropy has been studied extensively in the last few years. 
The basic goal has been to map its predicted I dependence in viable cosmological, 
large scale structure, and IC gas models. The strong motivation for this is the 
need to accurately calculate the power spectrum of the full anisotropy in order to 
make precise global parameter determinations from large stratospheric and satellite 
databases. In addition, mapping the S-Z anisotropy will yield direct information 
on the evolution of the cluster population. 

Results from many calculations of the predicted S-Z anisotropy are not always 
consistent even for the same cosmological and large scale structure parameters. 
This is simply due to the fact that the calculation involves a large number of input 
parameters in addition to the global cosmological parameters (e.g., the present 
cluster density, and parameters characterizing the evolutionary history of IC gas), 
and the sensitive dependence of the anisotropy on some of these. The anisotropy 
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can also be predicted - presumably more directly from simulations of the S-Z sky 
based largely on results from cluster X-ray surveys and the use of simple scaling 
relations (first implemented by Markevitch et al. 1992). 

The anticipated observational capabilities, to be achieved with long duration 
balloon-borne experiments and satellites, of detailed mapping of the small angular 
scale anisotropy, have motivated many recent works. Colafrancesco et al. (1997), 
and Kitayma et al. (1998), have calculated the S-Z cluster number counts in an array 
of open and flat cosmological and dark matter models, and Scljak et al. (2000) have 
recently carried out hydrodynamical simulations in order to generate S-Z maps and 
power spectra. Cooray et al. (2000) have, in particular, concluded that the planned 
multi-frequency survey with the Planck satellite should be able to distinguish be- 
tween the primary and S-Z anisotropies, and measure the latter with sufficient 
precision to determine its power spectrum and higher order correlations. 

The main characteristics of the predicted power spectrum of the induced S-Z 
anisotropy are shown in Figure 3. Plotted are the angular power spectrum func- 
tions, Ce,{l + l)/27r, vs. the multipole, £, for both the primary and S-Z induced 
anisotropies in the standard CDM model with ft = 1, based on the work of Sadeh 
and Rcphaeli (2001). Their treatment is an extension of the approach adopted by 
Colafrancesco et al. (1997), who used the Press & Schechter formulation for the 
calculation of the cluster density as function of mass and redshift, and normal- 
ized it at the current epoch by the observed X-ray luminosity function (see also 
Colafrancesco et al. 1994). IC gas was assumed to evolve in a simple manner con- 
sistent with the results of the EMSS survey carried out with the Einstein satellite. 
The primary anisotropy was calculated using the CMBFAST code of Seljak & Zal- 
darriaga (1996). The solid line shows the primary anisotropy which dominates over 
the S-Z anisotropy for I < 3000. The S-Z power is largely due to the thermal effect; 
this rises with £ and is maximal around £ ~ 1000. In this model, the S-Z power 
spectrum contributes a fraction of 5% (10%) to the primary anisotropy at £ ~ 1860 
{£ ~ 2360) . The relative magnitude of the S-Z power is higher in the A-dominated 
CDM model with A = 0.8; 5% and 10% contributions occur at £ ~ 1520 and 
£ ~ 1840, respectively. It can be concluded from this (and other studies) that the 
S-Z induced anisotropy has to be taken into account if the extraction of the cosmo- 
logical parameters from an analysis of measurements of the CMB power spectrum 
at I > 1300 is to be precise. 

5 Future Prospects 

Work on the S-Z effect has greatly advanced in the last few years, and the prospects 
are very good for major improvements in the near future. The full potential of the 
effect as an important cosmological probe will be realized when high quality spec- 
tral and spatial mapping of the effect in nearby clusters, z < 0.1, will be feasible. 
The main limitation on the accuracy of the cosmological parameters will continue 
to be due to systematic uncertainties, which can be most optimally controlled when 
sensitive S-Z and X-ray measurements are made of nearby clusters. The interfero- 
metric CBI array extended the capability of the BIMA and OVRO arrays to obtain 
sensitive maps of the effect in nearby clusters. The currently operational XMM and 
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Figure 3. Primary and S-Z power spectra in the flat CDM model (Sadeh & Rephaeli 2001). 
The solid (red) line shows the primary anisotropy as calculated using the CMBFAST 
computer code of Seljak & Zaldarriaga (1996) . The dashed (blue) line shows the thermal S- 
Z power spectrum, and the dotted-dashed (black) line is the contribution of the kinematic 
component. 



Chandra satellites provide state-of-the-art X-ray spectral and spatial capabilities. 
Use of the S-Z spectrum as a powerful diagnostic tool will soon be possible when 
sensitive multi-frequency bolometric array experiments, with spectral bands in the 
range 150 — 350 GHz, and very good spatial resolution (~ 3'), become operational. 
Both ground-based (e.g., the upgraded MITO experiment) and stratospheric (e.g., 
BOOST, OLIMPO) projects are underway. Measurements of a large number of 
clusters with these and the currently operational arrays will enable a much more 
precise determination of the cosmological parameters. In particular, it will be pos- 
sible to measure Hq with an overall uncertainty of just ~ 5%. 
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